基于多GPU的大型线性和非线性方程组的求解(英文)

在处理工程问题时,常常需要对线性或非线性方程组进行求解。对于实际应用中经常遇到的大型方程组进行求解则需要相当长的时间。使用图形处理器(GPU)代替传统的CPU,将多块GPU通过操作系统进行协调,并将PBi-CGstab方法和Inexact Newton方法进行适合多GPU并行的改造以此作为多GPU求解器的核心算法,加速求解大型线性和非线性方程组。本文的多GPU求解器在成倍扩展了单GPU求解器允许的计算规模的同时取得了令人满意的加速比。     

Modelling ignition probability with pairwise mixing-reaction model for flame particle tracking

Reduced order models for ignition analysis can offer insights into ignition processes and facilitate the combustor optimization. In this study, a Pairwise Mixing-Reaction (PMR) model is formulated to model the interaction between the flame particle and the surrounding cell mixture during Lagrangian flame particle tracking. Specifically, the model accounts for the two-way coupling of mass and energy between the flame particle and the surrounding shell layer by modelling the corresponding turbulent mixing, chemical reaction and evaporation process if present. The state of a flame particle, e.g., burnt, hot gas or extinguished, is determined based on particle temperature. This model can properly describe the ignition process with a spark kernel being initiated in a nonflammable region, which is of practical importance in certain turbine engines and has not been rigorously accounted for by the existing models based on the estimation of local Karlovitz number. The model is integrated into an ignition probability analysis platform and is demonstrated for a methane/air bluff-body flame with the flow and fuel/air mixing characteristics being extracted from a non-reacting simulation. The results show that for the spark location being at the extreme fuel-lean outer shear layer of the recirculation zone, PMR can yield ignition events with a significant number of active flame particles. The mechanisms for the survival of the initial flame particles and the entrainment of the survived flame particles into the recirculation zone are analyzed. The results also show that the ignition probability map from PMR agrees well with the experimental observation: a high ignition probability in the shear layer of the recirculation zone near the mean stoichiometric surface, and low ignition probabilities inside the recirculation zone and the top stagnation region of the recirculation zone. The parametric study shows that the predicted shape of the ignition progress factor and ignition probability is in general insensitive to the model parameters and the model is adequate for quantifying the regions with high ignition probabilities.     

ZL205A合金舱体低压工艺设计及数值仿真应用

针对某ZL205A合金舱体铸件壁薄难充型、壁厚差异大、易产生裂纹等特点,设计了低压底注式浇注工艺。采用立筒缝隙式浇注系统分配充型过程中合金的流量,局部厚大部位放置成形冷铁加快该部位的冷却速度,以控制铸造缺陷的产生。利用ProCast数值仿真软件对工艺进行仿真计算,成功预测了缺陷的产生位置,对工艺方案进行优化,包括倒角过渡、改变型砂使用等,并在优化方案的基础上浇注了舱体铸件,铸件质量良好,没有缩松、缩孔、裂纹等缺陷,符合相关技术要求。     

热塑性定型剂对干铺丝-RTM工艺及力学性能影响

为了研究热塑性定型剂的添加对干铺丝-树脂传递模塑（RTM）工艺的影响，采用PA-6作为热塑性定型剂，通过T-剥离确定定型剂的质量分数，结合试验和工艺仿真，研究了定型剂对预制体面内渗透率的影响，并对其充模过程进行工艺仿真和分析，在此基础上研究了定型剂对环氧复合材料制品层间性能及韧性的影响。渗透率试验和仿真结果表明:加入定型剂后，纤维预制体的沿纤维和垂直纤维的面内渗透率分别降低了55.5%和52.0%；工艺仿真结果与试验结果一致，不同渗透长度下的注胶时间偏差均小于10%。力学性能试验结果表明:加入热塑性定型剂可显著提升环氧树脂复合材料的韧性，Ⅰ型和Ⅱ型层间断裂韧性分别提高了89.3%和54.1%，冲击韧性提高了18.7%，复合材料的层间剪切强度降低了23.6%。      

Sliding mode control design for oblique wing aircraft in wing skewing process

When the wing of Oblique Wing Aircraft (OWA) is skewed, the center of gravity, inertia and aerodynamic characteristics of the aircraft all significantly change, causing an undesirable flight dynamic response, affecting the flying qualities, and even endangering the flight safety. In this study, the dynamic response of an OWA in the wing skewing process is simulated, showing that the three-axis movements of the OWA are highly coupled and present nonlinear characteristics during the wing skewing. As the roll control efficiency of the aileron decreases due to the shortened control arm in an oblique configuration, the all-moving horizontal tail is used for additional roll and the control allocation is performed based on minimum control energy. Given the properties of pitch-roll-yaw coupling and control input and state coupling, and the difficulty of establishing an accurate aerodynamic model in the wing skewing process due to unsteady aerodynamic force, a multi-loop sliding mode controller is formulated by the time-scale separation method. The closed-loop simulation results show that the asymmetric aerodynamics can be balanced and that the velocity and altitude of the aircraft maintain stable, which means that a smooth transition is obtained during the OWA’s wing skewing.     

基于产品结构的航天器AIT过程多维数据管理系统

文章为解决航天器AIT过程中数据集成度低、数据应用不佳等问题,从数据模型构建、数据集成等方向出发,分析了AIT过程数据在产品结构、时间周期、阶段型号3个维度的演变过程,构建了多维数据模型,并开发基于产品结构的多维数据管理系统。该系统以产品结构为核心,对AIT过程数据进行重新组织,定义不同系统间、不同维度数据之间的关联关系,实现了AIT过程全数据的结构化整合,在此基础上生成了AIT过程数据包。     

航空涡扇发动机起动热悬挂故障诊断

为了完善某型航空涡扇发动机起动性能设计,通过分析发动机的起动控制逻辑和热悬挂控制逻辑,确定了起动热悬挂故障的原因。起动热悬挂分为排气温度 T6上升率过高故障和压气机出口压力 P31上升率过低故障 2种模式,故障原因分别为初始供油量偏多和起动过程 P31波动,针对故障原因制定了解决措施:初始点火供油量设定值下调,提高 P31传感器小压力段的测量精度,并优化 P31斜率异常故障判故逻辑。上述解决措施不仅解决了起动热悬挂问题,更进一步增强了产品设计与生产的兼容性,提高了发动机的起动性能和起动成功率。     

三维五向编织复合材料的冲击压缩特性及破坏机制

利用分离式霍普金森压杆(split Hopkinson pressure bar,SHPB)实验系统和高速摄像技术对三维五向碳/环氧编织复合材料的动态压缩特性进行研究。通过对编织角为22.3°的试样分别进行沿纵向和横向方向的冲击压缩实验,得到材料在200~1200 s^-1应变率范围内的应力应变曲线,并结合高速摄像记录的动态压缩过程,对不同应变率下材料在高速变形下的渐进破坏规律进行分析。同时,综合试样的宏观破坏特征和微观断口形貌特征,进一步分析材料的破坏模式及破坏机理。结果表明:随着应变率的增加,材料在纵向和横向均具有一定的应变率强化效应,在横向方向的应变率强化效应更为显著;不同加载方向下材料的渐进破坏过程、应力应变曲线特征以及破坏方式均具有明显差异,且随着应变率的增加而发生改变。     

流量调节器管路系统自激振荡特性研究

流量调节器管路系统在小流量大压降工况下会出现低频自激振荡现象。为了深入认识自激振荡产生机理,结合某稳流型流量调节器及管路系统,基于流量调节器弹簧振子动力学模型开展数值仿真研究。数值仿真得出自激振荡频率为94Hz,与发动机试验结果一致。分析了流量调节器结构参数对系统稳定性的影响作用,三角形滑阀节流口能够抑制管路系统自激振荡。自激振荡产生机理是液动力随滑阀节流口型面振荡,并对管路系统形成正反馈作用,当流量调节器综合刚度系数<0时,管路系统就失稳产生振荡。某流量调节器的负载特性试验表明,随着流量调节器压降升高,管路系统稳定性变差。仿真获得的幅频特性,稳定边界与试验结果一致。     

流量调节器管路系统自激振荡特性研究（液体火箭推进专刊）

流量调节器管路系统在小流量大压降工况下会出现低频自激振荡现象,为了深入认识自激振荡产生机理,结合某型自稳流型流量调节器及管路系统,基于流量调节器弹簧振子动力学模型开展数值仿真研究。数值仿真得出自激振荡频率为94 Hz,与发动机试验结果一致。自激振荡产生机理是流量调节器的液动力受滑阀窗口型面变化率影响对系统形成正反馈作用,流量调节器综合刚度小于零时系统失稳。分析了流量调节器结构参数对系统稳定性的影响作用,三角形滑阀节流口流通面积能够抑制管路系统自激振荡。结合某型流量调节器负载特性试验系统进行验证,得出随着流量调节器压降升高,管路系统稳定性变差,获得的幅频特性,稳定边界与试验结果一致。